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ABSTRACT 

To ascertain whether magnetic dynamos operate in rocky exoplanets more massive or hotter than 
the Earth, we developed a parametric model of a differentiated rocky planet and its thermal evolution. 
Our model reproduces the established properties of Earth's interior and magnetic field at the present 
time. When applied to Venus, assuming that planet lacks plate tectonics and has a dehydrated mantle 
with an elevated viscosity, the model shows that the dynamo shuts down or never operated. Our model 
predicts that at a fixed planet mass, dynamo history is sensitive to core size, but not to the initial 
inventory of long-lived, heat-producing radionuclides. It predicts that rocky planets larger than 2.5 
Earth masses will not develop inner cores because the temperature-pressure slope of the iron solidus 
becomes flatter than that of the core adiabat. Instead, iron "snow" will condense near or at the top 
of these cores, and the net transfer of latent heat upwards will suppress convection and a dynamo. 
More massive planets can have anemic dynamos due to core cooling, but only if they have mobile lids 
(plate tectonics). The lifetime of these dynamos is shorter with increasing planet mass but longer with 
higher surface temperature. Massive Venus-like planets with stagnant lids and more viscous mantles 
will lack dynamos altogether. We identify two alternative sources of magnetic fields on rocky planets: 
eddy currents induced in the hot or molten upper layers of planets on very short period orbits, and 
dynamos in the ionic conducting layers of "ocean" planets with ^-^10% mass in an upper mantle of 
water (ice). 
Subject headings: planets and satellites; interiors, magnetic fields, tectonics, physical evolution 
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1. INTRODUCTION 

Earth distinguishes itself in the inner Solar System 
with a magnetic dipole generated by a dynamo in its 
iron core. Mars had a magnetic field of comparable 
intensity that collap sed early in the plan et's history 
(I Arkani-Hamedl[2004[ ) and never re-appeared (jLillis et al.l 
I2008D. Mercury ha s a weak (1% of E arth) global field 
(jStanlev et aIll2005l:[Christensenll2006[ ). Venus currently 
lacks a dynamo, and we have n o evidence for or against 
one in its past (jNimmol 120021 ). At least one massive 
sateUite in the outer So lar System (Ganymed e) appears 
to have a core dynamo ([Schubert et al.|[l996D . 

Planets as small as twice the mass of the Earth 
have been discovered on shor t-period orbits around 
other stars (jMavor et al.l 12009^. A fraction of these 
transit the parent star, allowing mean densities to be 
determined. The smallest transiting exoplanet has a 
mass of 4.8 ± 0.8 M^ and a mean density of 5.6 ± 1.3 
g cm~^, consiste n t with an Earth-l ike composition 
(|Leger et al.) 120091: IQueloz et all 120091 ). The second 
smallest (6.6 ± 0.9 M^) has a mean density only 
one-third that of Earth, indicating it has a thick volatile 
envelope (|Charbonneau et al.l 120091 ). These two planets 
may rotate synchronously and have effective emitting 
temperatures of around 2600 K and 500 K, respectively 
(that of Earth is 255 K). The Kepler spacecraft mission 
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(jKoch et al.l 120101: IBorucki et all 120101) is expected to 
find many such "hot" ma ssive rocky planets ([Selsis et al.l 
[20071: iGaidos et al.l[2007l) . 

Do these planets have magnetodynamos? While there 
have been many theoretical studies of the dynamos in 
the Earth and smaller So lar System bodies, as well as 
in th e ice and gas giants ([Breuer et al.l 120091 : lStevenson[ 
120091) ■ rocky planets more massive than the Earth have 
not been considered, although simple dynamo scaling 
laws have been applied to ga s giants on short-period 
orbits (Greif5meier et al] I2004D . Magnetic fields may 
protect planetary atmospheres against erosion by stellar 
wind s and coronal mass ejections (iGreiBineier et al.l 
120041: iLundkin et al.ll2007l: iDehant et al.ll2007[ ). Interac- 
tion between a giant planet's magnetic field and that of 
the host star has be en inferred from corotational chro- 
mosp heric activity ([Shkolnik et al.l I2008t IWalker et al.l 
I2008D and may produce detect able radio emission 
( Lecavelier Pes Etangs et al.l 120091) . Ohmic dissipation 
of the kinetic energy of winds in partially ionized 
atmospheres has been proposed as an explanation for 
the "inflated" radii of some short-period gas giants 
([Batvgin fc Stev enson 2010). 

For a planet to have a dynamo, it must contain a fluid 
layer that is sufficiently electrically conducting for the 
magnetic Reynolds number Re™ = VL/\ > 40, where 
V and L are a characteristic fluid velocity and length, 
respectively, A = \/{^Qa) is the magnetic diffusivity. 
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/io = 47rxlO^^NA^'^, and a is the electrical conductiv- 
ity. In Earth-like planets the conducting flui d is liquid 
iron or an alloy thereof (A ~ 2 m^ sec~^ iStevensoiil 
([2009)). The possibility that the oxidation state of 
planetary material might prev ent iron co re formation 
was explored by Elkins-Tanton fc Seageil p008) . but 
we assume that such cores do form. Metallization of 
silicates and miscib ility with Fe is not expected until a 
pressure of 20 TPa (jOganov et al .1120051) . In the cores of 
Earth-size planets (L ~ 3 x 10^ m), convective motions 
as small as 10~^ m sec~^ are sufficient for Re™ > 40. 
The temporal variation of Earth's magnetic field, such as 
westward drift, imply velocities of 10~* — 10""^ m sec~^ 
at the top of the core. 

A dynamo also requires a source of energy to maintain 
convective motions against internal ohmic dissipation. 
Core convection can be driven by the removal of heat 
from the core faster than it can be transported by 
conduction, by the release of latent heat during phase 
transitions (i.e., solidification of iron), or by the forma- 
tion of buoyant fluid via the expulsion of light elements 
during iron solidiflcation. The last two are related to 
the first because the latent heat of solidification must 
be removed. Thus dynamo activity is ultimately related 
to the transport of heat and the temperature contrast 
across the core-mantle boundary (CMB). That contrast 
can be maintained either by efficient cooling of the lower 
mantle (e.g., by descending slabs of former lithosphere 
during plate tectonics) or by the heating of the core by 
radioactive nuclides. The absence of a dynamo in Venus 
has been ascribed to the lack of plate tectonics, which 
slows cooling of the mantle, as well as the core (,Nimmoi 
[2003). 

In dynamo theory, the condition for convection is ex- 
pressed as a requirement on entropy production, rather 
than energy production, because the latter is unaffected 
by dynamo dissipation when mechanical energy is con- 
verted to heat. The balance between entropy sinks and 
sources in the core gives the maximum dissipation $ by 
the dynamo. If $ > then a dynamo is permitted. We 
calculate $ using parameterized models of the interior 
structure and thermal evolution of a planet consisting 
of an Fe core and a silicate mantle. We account for, 
but do not vary, the presence of light elements in the 
core. We estimate the magnetic Reynolds number and 
average surface field using previously established scaling 
laws. We address whether dynamos operate in rocky 
planets with different masses, surface temperatures, and 
modes of mantle convection. Figure [T] is a schematic 
of our simplified planet. In the equations that follow, 
we employ the following subscripts: 1 (ambient), p 
(planet or surface), m (mantle), c (core or core- mantle 
boundary), i (inner core boundary), and (center). 

2. INSIGHT FROM SCALING LAWS 

The manetostrophic approximation assumes that 
Lore ntz forces are b alanced by Coriolis forces in the 
core (lStevensonll2003D . This produces a scaling relation- 
ship between magnetic field strength (or dipole moment) 
and core mass, radius, and rotation rate. Alternatively, 
ohmic dissipation, and hence the strength of the dynamo. 



are s et by the available power (jChristensen fc Aubertj 
120061) . Analysis of numerous numerical dynamo simu- 
lations shows that, in this regime, the RMS strength of 
the field in the core, Bc^ can be related to a Rayleigh 
number based on a mass anomaly flo w (mass deficit 
advected per unit time) (iChristensen fc Aubertl I2Q06I : 
lAubert et a Ill2009l:ICh ristensenl l2Q 09D. and thrmig h that 



to the convective power in the core (jBuffett et allll996D . 
In the case where the magnetic Prandtl number is much 
less than one and essentially all the convective power is 
lost to ohmic dissipation, that relationship can be ex- 
pressed in terms of a dimensionless convective power p 
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where ai and 6i are dimensionless parameters, p is the 
mean density of the core, n is the angular rotation rate of 
the planet, Ro and Ri are the outer and inner boundaries 
of the convecting zone, and p = (j)T / \p?{Rc — Rif\. 
(pT is the available convective power per unit mass, 
where cj) — ^ /M^ is the entropy available per unit mass 
and time, and T is an effective dissipation temperature. 
lAubert et al.l (|2009[ ) find oi ~ 1.65 and 6i to be almost 
exactly 1/3. The latter eliminates the dependence on Vt: 



/Jm^ [(t)T{Rc - Ri 



I 1/3 
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Likewise, the magnetic Reynolds number can be found by 
a relationship between the Rossby number and p, giving: 
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where again a2 and 62 are dimensionless. 62 is not nec- 
essarily 1/3 and t his introduces a sligh t dependence on 
VL. For example, lAubert et all (|2009t ) find 02 « 1.3 
and 62 ~ 0.42. For the purpose of a simple order-of- 
magnitude estimate of Re™, we ignore these complica- 
tions, set 62 = 1/3, producing a new formulation: 
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Re — Ri 
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(4) 



For plausible terrestrial values (0 ^ 100 MW K ^, 
f - 5000 K; iLabrossd ([200I ). Re™ ~ a^lOl So as 
long 02 ^ 1 its exact value will not be critical. Com- 
bining Equations [2] and [4] for terrestrial values gives 
Re™ ~ 16(i3c/l/^T). Thus if the predicted field in the 
core is at least a few /iT, i.e. the surface field (attenuated 
at least by the cubic distance law of a dipole) is > 1/iT, 
then Rcm > 40. This criterion is only very weakly de- 
pendent on planet mass. 

The strength of the field at the surface of the planet 
and beyond will be sensitive to the field's topology, es- 
pecially the fraction in the dipole component, and will 
depend on fi. More rapidly rotating planets (lower 
Rossby number) will have a stronger dipole component 
(jChristensen fc AubertI I2006D . In the case of a pure 
dipole, the average surface field will be 



Bp ~ Bc{Rc/RpY 



(5) 



We estimate the Bp and Re™ using Equations [2l |4l and 
[SI We conservatively take T to be the temperature at 
the top of the core. We assume that the field is dipole- 
dominated like the Earth. 
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3. MODEL 

3.1. Interior structure 

A planet is modeled as a homogeneous, fully convect- 
ing Mg/Fe-silicate mantle surrounding a liquid/solid Fe 
core. The pressure-density profiles are calculated self- 
consistently with the (time-varying) size of the inner 
core. Third-order Birch-Murnagham (BM) equations of 
state (EOS) are used for each component. The pressure 



P = -K, (x^ 

2 ^ 



-{4-K[){l-x' 



(6) 



where x — {p/ piY^^ is the isotropic strain. Values of the 
ambient density pi, compressibility Ki and its pressure 
derivative K[ for the relevant planetary materials are 
given in Table 1. Electron degeneracy (Thomas- Fermi- 
Dirac, or TFD) pressure will be important in the metal 
cores of massive Earth-like planets. We calculated the 
TF D c ontribution to the pressure using the formulation 
of IZapolskv & Salpctcr (1969) and found that the 
EOS for liquid and e-phase solid iron (Fe(l) and Fe(e)) 
intersect at about 3.3 TPa. Pressures beyond 3 TPa 
are not relevant to the range of masses considered here. 
The density in the liquid part of the core is adjusted by 
a fixed fr a.ction 5pl p to account for the presence of light 
elements (|Li fc Feil 120031 ) . 

Our interior model reproduces the central pressure 
(364 GPa), CMB pressure (1 36 GPa), and core size 
Rr ( 3480 km) of the Earth (jDziewonski fc AndersonI 
Il981h . but only if the zero-pressure compressibility for 
Fe(l) is 40% larger than th e 106 GPa value reported by 
[Anderson fc AhrensI (|1994( ) based on shock experiments. 
With this adjustment, the predicted radius of the Earth 
is 25 km (0.4%) less than the actual value, but this is 
not surprising because we do not include low density 
crustal rocks or an ocean. Our dynamo predictions are 
much more sensitive to the properties of the core than 
those of the entire planet. Predicted planet radius, core 
radius, pressure at the CMB, and central pressure are 
plotted versus total mass in Figure [2l Pressures increase 
as the core solidifies and contracts. 



3.2. Temperature and density profiles 
The adiabatic temperature profile in the outer convect- 



ing core is approximately 

T(r)«T,exp[(i?2_ 



')IDI 



(7) 



where the thermal length scale evaluated at the planet 
center is, 



Dn 



OC-n 



(8) 
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(jLabrossd I2003D , Cp and a are the heat capacity and 
thermal expansivity of Fe, respectively, and G is the 
gravitational constant. The temperature profile deviates 
from Equation[7]to the extent that D, and specifically a, 
vary with depth in the outer core, a must therefore be 
evaluated appropriately when calculating terms that are 
sensitive to I?, i.e. the heat and entropy transported by 
thermal conduction. Hereafter D (without a subscript) 



is Dq. 
The density in the core follows 

p{r)^p,exp[{Rl^r^)/L'] 
where the density scale length is 



L = 



9Ki 



2ttGpiPo 



In^ 
Pi 



1 



(9) 



(10) 



For the Earth, D -- 6400 km and L -- 7400 km and these 
values are only weakly dependent on planet size. We 
expect the core radius i n. planets of Earth- like compo- 
sition to scale as M^-^s (jSotin et al.l 120071 : iSeager et all 
l2007f l. The nth order term in either r/D or r/L will be 

of order Mp /2". For planets of a few Earth masses, 
we must retain terms in expansions with r/D or r/L 
through fifth order. 

If a solid inner core is present, the temperature at the 
inner-outer core boundary is the intercept between the 
adiabat and the iron melting point r. Adopting a Lin- 
demann law for r. 



51ogT 
dlogp 



= 2 
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(11) 



where 7 is the Griineisen parameter. The application 
of the Lindemann law to the deep int eriors of planets 
issometimes viewed as speculative (e.g.. lWolf fc Jeanlod 
()1984() ). The assumptions upon which it is derived are 
not strictly valid for polyatomic systems, and the law is 
known to sometimes fail for minerals with complex in- 
teratomic forces and structures. However, at high pres- 
sures and for metals it becomes a better approximation, 
and hence is widely used to interpret shock data (e.g., 
[Anderson fc AhrensI ()1996D V Importantly, it appears to 
provide a good description of expe rimental da ta and first 
princip le calculations for Fe fe.g.. iHenilev fc Maa (12001) : 
IHuseri (|2007D ) and other metals (e.g. IDai et al.l (|2002[) 1. 



The radial dependence of the solidus becomes 



T{r) = To exp 



-2 1- 



!7 j \d) 



Si J 



(12) 



where we have used L^/D^ = 7 (|Labrossell2003[) . Thus, 
the temperature at the top of the core is uniquely set by 
the ratio of the inner to outer core radius SR = Ri/Rc- 



Tc = tq exp 



1 + 3fiM 1 



2 

3^ 



RP 



(13) 



Tc must decrease as the inner core grows. If the entire 
core is liquid, then Tc is the independent variable. 
For greater accuracy, we use our interior model to 
directly calculate the density at the ICB, the melting 
temperature using Eguation llll and Tc using Equation[71 

The temperature profile of the inner core will de - 
pend on whether it, too, is convecting (jBuffetti [2009D . 
Such convection may be sufficiently torpid so as not to 
contribute to the dynamo, but as long as overturn is 
rapid compared to the timescale of inner core growth 
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(10^ — 10^" yr), the temperature profile will be an adi- 
abat. This assumption maximizes the transport of sen- 
sible heat into the outer core, and hence minimizes the 
rate of inner core growth. 

For a fully convecting core, the heat conducted along 
the the adiabat at the CMB is (neglecting compressibility 
effects) 

AirR'^kacgcTc 



Q 



K 



(14) 



where k is the thermal conductivity and g is the local 
gravitational acceleration. If Qc, the heat flow from 
the core to the mantle across the CMB, is less than 
this, the temperature profile will be sub-adiabatic to the 
degree required to equalize the heat flows. Whether this 
layer is completely stratified will depend on non-thermal 
(i.e., light element) buoyancy forces. In the absence 
of a growing inner core, however, such a layer will be 
thermally stratified. Geomagnetic variation limits the 
thickness of any upper stratified layer in the Earth's 
core to < 100 km (Gubbins 2007i) , if one exists at all 
(jStanlev fc Mohammadill2008D . 

If the stratified layer is thin and makes a negligible con- 
tribution to the sensible heat fiow, the heat fiow through 
it is constant and its temperature profile is simple. The 
location 5R, and temperature T* of the boundary between 
the convective and conducting regions is found by simul- 
taneously solving for the heat flow and temperature: 



m^ 



SirkDrT^ 



and 



1 c — -^ * 



Qc{Tc 
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and using Equation [71 i.e. 



T* = Ti exp 



{^'~^l 



- 1 



(15) 



(16) 



(17) 



where the dependence of Qc on Tc must be accounted 
for ( Ji3.4p . and we evaluate the thermal scale length at 
the CMB. In the event that stratification does occur, 
we compute the entropy terms for the convecting part 
of the core only. This is done by re-scaling the radius 
and mass of the core in the entropy equations to the 
size and mass of the convecting zone. T^ becomes the 
temperature T* at the top of the convecting zone, and we 
assume that since the adiabat is self-similar, such that 
dT^,/T^, « dTc/Tc (effectively ignoring small differences 
due to changes in the thickness of the convective zone 
itself). 



3.3. Equation of entropy production 

The condition for maintenance of convective motions 
for a dynamo can expressed as a bal ance between so urces 
and sinks of entropy in the core, e.g lLabrossd ()2007D and 
iNinmiQ (2009): 



^ + Ek ^Er + Es + Eg + El, 



(18) 



where Ek is the entropy sink due to conduction along 
the adiabat, Eji the production of entropy by internal 
heating; Es that due to core cooling (loss of sensible 
heat); Eg the entropy production from buoyancy gener- 
ation; and El the entropy production due to latent heat 
generated by crystallization of the core. The decay of 
radioisotopes is a possible internal h eat source. ^ °Fe will 
decay completely (ii/2 = 2.6 Myr) ()Ruge]| l2009t ) before 
planets finish accreting, but '"'K (^ 1/2 = 1-25 Gy r) may 
be a heat source in Earth's core (jNimmo et al.l 12004 ). 
Values of the partition coefficient of K between silicate 
and Fe/Ni alloy melts have been discrepant between dif- 
ferent hi gh pres sure experiments, possibly for technical 
reasons (|Li fc F ei 2003), but new results limit the abun- 
dance of K to a few tens of ppm, wi th no evidence for a. n 
increase in solubility with pressure (jCorgne et al.ll200"7[ l. 
In the absence of radiogenic heating, 

<^ = El+Es+Eg-Ek. (19) 

We calculate $, divide by the outer core mass to obtain 
(/), and use Equations[21[SJ and|l]to estimate the magnetic 
field intensity and magnetic Reynolds number. 

Expressions for entropy production (dissipation) have 
been deri ved and disc ussed exte nsively in the li tera- 
ture, e .g. Listed (|2003t ). Following Eabrossd (|2007| ) and 
iNimmo (2009) we derive expansions in terms of the size 
of the inner core. The term from the loss of sensible heat 
comes form the thermodynamic identity dS = CpdT/T 
under constant pressure and no heats of reaction. Thus 

^-i)po,fdF, ,20) 



Es = - 



where the integration is over the volume of the entire 
core. In the absence of an inner core, and to fifth order 
in Re/ D and Rc/L, 
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-CpMc 
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7\D 
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dT^ 
dt 



(21) 

where Mc is the total mass of the core. We assume 
that subsolidus convection continues the adiabatic 
temperature profile into the inner core ([Buffettl |2009() 
and ignore the effect of small changes in the temper- 
ature and pressure length scales due to differences in 
compressibility and density between the liquid and solid 
(e) phases. This assumption maximizes the contribution 
of the inner core to the entropy and energy terms. In 
the case of a mostly solid core, a detailed treatment of 
the heat flow from the inner core would be warranted. 
However, we shall show that in planets much larger than 
the Earth, inner cores remain relatively small, or do not 
grow at all. 

The cntropic term for the latent heat due to core 
growth is; 



El = 47ri?,^p,A5 ( ^ _ 1 ) ^ 



Tr 



dt 



(22) 



where AS" is the specific entropy of fusion. Assuming 
that the core grows outwards. 



dRi 

~dF 
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2Ri (A - 1) Tc dt ' 



(23) 
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where A is the ratio of the temperature-pressure slopes 
of the sohdus to the adiabat, evaluated at Ri. Then, 
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3 McAS Pi 
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SR d - 3ff2 



l-3f{2 



Tc dt ' 

■(24) 



where the ratio of the density at the top of the inner 
core to the mean density is obtained directly from the 
interior model ( ^3.1[) . 



The entropic term for the release of buoyant fluid (den- 
ity deficit Ap 
(|Labrossdl200i 



sity deficit Ap) during the crystallization of the core is 

]3); 
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dt 
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where Ap is the density change in the fluid due to the 
presence of light elements, and F and F' are the dimen- 
sionless functions 
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c.f iNimniol (|2OO90 . To derive this, we ignore work done 
by expansion of the core (Li ster 2003 ), and assume that 
the energy is dissipated at temperature Tj,. 

The entropy production by the conduction of heat 
along the adiabat is 



E 



K 



Airr^k i Tf^-^ \ dr 



(|Listeiil2003l ). This is evaluated to be 



Ek = 



12kMc 



l-^" 



(27) 



(28) 



To more accurately reflect the temperature profile, to 
which conduction is sensitive, we have evaluated the 
temperature scale length at the CMB, rather than the 
core center. The other terms in the entropy equation 
depend on £> to a lesser degree and are evaluated over 
the entire core, not just the CMB. 



3.4. Core cooling rate 
We solve for the core cooling rate dTc/dt by balancing 



the heat fiow: 



Qs + Ql+Qg = Qc, 



(29) 



where the terms correspond to sensible heat, latent heat, 
gravitational energy, and heat transported through the 



CMB, respectively. The left-hand terms are all propor- 
tional to dTc/dt. The sensible heat term is 

Qs^ j P^P^dV. (30) 

Assuming the adiabat continues into the inner core, this 



Qs 



CpMc 



4 
35 



The latent heat term is 

3 MrAS^ f D_ 
R'r 



175 \DL 

(31) 
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2 A- 1 
while the gravitational term is. 



exp[iRjD)\l-di^)] 



dT^ 
dt 



SirGMcApD^ p, F dT^ 
A- 1 'pYc'~dt 



(32) 
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We compare Qc to the heat Qk that can be carried by 
conduction along the core adiabat: 



!K 



%TTRlkTc 



(34) 



Qc = AirRchnNuc (Tc ~Ti), 



(35) 



where T/ is the temperature of the lower mantle and km 
is the thermal conductivity of the mantle. The Nusselt 
number is 

Nuc ~ (Rac/Ra*)^ (36) 

where 5 « 0.3 (|Schubert et al.ll2001h . For a fluid heated 
from below, 

pga {Tc - Ti) d^ 



Rac = 



KTjc 



where the gravitational acceleration g and viscosity rj are 
evaluated at the CMB, and d = Rp — Re- In the case of a 
partially stratified layer, we substitute Tt for Tc, where Tt 
is the temperature at the top of the stratified, conducting 
zone. 

The heat transport by subsolidus convection in the 
mantle is; 



Q. 



ATrRlNumk (T^ - Tp) 



Rv — Re 



(38) 



dTc 
dt 



If Qc > Qk then the entire core has thermally driven 
convection. If Qc < Qk then cooling is insufficient to 
maintain convection in the entire core and a stratified, 
conducting layer develops ( ii3.2[) . In our calculation of 
the rate of sensible heat loss, we have neglected the 
small corrections due to departures of the temperature 
profile from the purely adiabatic case. 



3.5. Heat transport in the mantle 

We assume that the mantle boundary layer at the 
CMB develops independently of the flow in the rest of 
the mantle and its thickness is controlled by a local 
Rayleigh numb er Ra, as oppo sed to the whole-mantle 
Ra (Nimmo fc Stevensonll2000h . Heat is removed from 
the core according to; 
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where Tp is the surface temperature. For a uniform in- 
ternal heat source, 



K,a,T| — 



pgHadP 



(39) 



where H is the total specific heat generation (per unit 
mass), rim is evaluated at the mantle reference tempera- 
ture, taken to be halfway along the adiabat. We make the 
approximation that heat from the CMB is mixed rapidly 
and uniformly into the mantle (e.g., by plumes) and thus 
H — He + Hr, where He = Qc/Mm, where Mm is the 
mass of the mantle, and the radiogenic heat production 
Hr is given by; 



lir — ^i^i^^i^ 



't/ti 



(40) 



Ci is the initial concentration of the ith radioactive iso- 
tope, hi is the specific heat production, and U is the 
mean life. The concentrations adopted for the long-lived 
radioisotopes '^''K, ^'^^Th, ^^^U, and -^■^^ U are for the case 
of an "undepleted" terrestrial mantl e (|Ringwoodlll99lD 
given in Table 2 of iKite et all ()2009[ ). The contribution 
due to mantle cooling is; 



Qi 



,M„ 



' dt 



(41) 



where Cm is the is the specific heat capacity of the mantle. 
We adopt the viscosity law for a stress exponent 
n = 1 material, appropriate for the plate-tectonic regime 
(jNimmo fc Stev enson 2000): 

Vrn = ■>]* exp {bTrn/T) , (42) 

where fe w 17 (iKarato et al.ll200ll ). t^ is the peridotite 
solidus, and 77* is a reference viscosity. In Equation 
[37l we evaluate ij at the "film" temperature, which is 
inter mediate that of th e lower mantle and top of the 
core ()Mangaet al.ll200lD . 

Convection in the silicate mantle may involve a 
mobile lid (i.e., plate tectonics) or a stagnant lid, and 
in our Solar System we have Earth-mass examples 
of each (Earth and Venus). Whether one mode or 
another will be more likely on larger planets will de- 
pend on the relative magnitudes of lithosphere stresses 
(iValencia et al.ll2007f). the ayailab litv of water to weaken 
faults iNimmo fc McK enzigll998f ). surface temperatures 
(JLenardic et al. 2008) . and the forma tion of buoyant, 
unsubductable crust (JKite et al.l 120091 ). In our model, 
the mode of mantle convection controls (a) the tempera- 
ture boundary condition at the mantle side of the CMB; 
and (b) the Nussult-Rayleigh number relationship in the 
mantle. 

In the case of a mobile lid, we assume that subducted 
slabs reach the lower mantle, have a temperature equal 
to the mean mantle temperature Tm, and directly cool 
the core. Thus for Equation [35] we set Ti — r,„. We the 
use the same Nu-Ra relation as in Equation 1361 In the 
stagnant lid case, the lower mantle temperature is set to 



Ti=9 exp 



" ag 



dr 



Ra ^P 



(43) 



and the relationship between mean and potential tem- 
peratures is taken to be 



= Tmcxp 



ag, 
— dr 



flc 



(44) 



We adopt the relationship between Nusselt number, 
Rayleigh number, and viscosity contrast N across the 
conducting lid (Frank-Kam enetskii parameter) given by 
iSolomatov fc Moresil (|2000( ): 



Nu = aN-^^+^^Ra^' 



(45) 



where a ~ 1, and boundary layer stability analysis shows 
that for a materia l with stress exponent n, /3 — n/(n + 2) 
(ISolomatovll 19951) . We use /3 = 0.6. 



4. RESULTS 

4.1. Parameter values and a nominal Earth-mass case 

Our numerical experiments revealed that the predicted 
history of the dynamo in a planet with an Earth-like 
interior is sensitive to: (1) mantle viscosity; (2) the 
thermal expansivity of iron (related to its compressibil- 
ity); (3) the thermal conductivity of Fe; (4) initial core 
temperature; and (5) the heat flow across the CMB. 
The first is constrained by the requirement to reproduce 
the size of and pressures in the Earth's core f ^3.1l) . We 
parameterize the last by using different values of the 
critic al Rayleigh number R a ^ . Alt h hough Ra„ is of order 
10^ JNimmo fc Stevensm] 120001 : iTurcotte fc Schubert! 
|2002| ) , we consider a range of values that encapsulate the 
complex uncertainties of the thermochemical boundary 
layer at the CMB. 

Mantle viscosity: Our viscosity law requires a solidus 
and a reference viscosity. W e fit a Lindemann melting 
law to the solidus data of iZerr et aD ()1998[ ) for the 
lower mantle, using the Mg-perovskite equation of 
state to relate P to p: We found a best-fit Griineisen 
parameter of 1.45. The reference viscosity of a hydrated 
Earth-like mantle was adjusted to achieve a present 
potential temperature 6 of 1750 K. This is the present 
temperature beneath mid-ocean ridge basalts (1623 K) 
(jHerzberg et al.l[201Cl( ) to which has been added 130 K 
associated with the olivine-spinel transition in the 
deep mantle (latent heat of reaction = 160 kJ kg~^) 
(|lto et al.l [19710. The derived reference viscosity, 
7y* = 6 X 10^^ Pa s, is somewhat above inferred values for 
the terr estrial asthenosphere where partial melting may 
occur (Solomatov & Moresi' '2000"; 'Dixon et al.' I2004J: 
F^eldskaar & Cathles 2007; BiUs ct al. 2007; James et al.l 
2009D . The viscosity of a dehydrated mantle is expected 



where the 6 is the potential temperature of the mantle. 



to be significantly higher. The effect of water on vis- 
cosity depends on the dominant creep mechanism, and 
the partitioning of water amongst the mineral phases 
that are present. Hydration of olivine in the upper 
mantle could l ower the viscosity by a factor of ^500 
(iHirth fc Kohlst edt 1996), but the degree of hydration, 
and concomitent decrease in visco sity, could be much less 
in th e lower perovskite mantle (jBolfan-Casanova et al.l 
I2OOOD . We selected a value based on the requirement 
that a simulation of Venus (0.815 M^ and Tj, = 785 K) 
predicts no dynamo at the present time. [The slow 



MAGNETODYNAMOS IN ROCKY EXOPLANETS 



rotation of Venus cannot be solely res ponsible for 
the absent dynamo (jZhang fc ZhandflQQSD .] If Venus 
has had a stagnant lid over its entire history, then 
the required viscosity enhancement for the Venusian 
dynamo to cease by 4.5 Gyr is a factor of 6. If its 
lid was primarily mobile, then the required viscosity 
enhanceent is a factor of 200. We conservatively adopt 
a factor of 10, recognizing that higher values are possible. 

Fe thermal conductivity: Thermal conduction com- 
petes with advection in removing heat from the core. 
Higher thermal conductivity will suppress convection, 
and thus the dynamo, and will stratify some or all 
of the liquid core. The thermal conductivity of iron 
is expected to depend on pressure, temperature, and 
the presence o f light elements, but not on temperature 
(jBermanl 119761 ) . The value often used in core models is 
50-60 W m-i K (|Labrossell2003f) . Thermal conductivity 
can be determined from measurements of electrical con- 
ductivity using the Wiedemann-Franz-Lorenz (WFL) 
law. Assuming the WFL law applies at high pres- 
sure, electrical con ductivity measureme nt at 140 GPa 
(jKeelerfeMitchelll [19691 and 200 GPa (|Bi et al.ll200^ 
imply k ^^100 W m"^ K-\ ci ., lAndersonI (fl998h . In 
contrast. iStacev fc Lopeii (|2007f ) argue for a value of 28- 
29 W m~^ K~^ based on the expectation that scattering 
of electrons by 3d orbitals (which contributes to resistiv- 
ity) will increase faster with pressure than the densities 
of 3d and 4s states (which contribute to con ductance). 
The e lectrical conductivity measurements of iMatassovl 
(|1977| ) show that light elements can sig nificantly reduce 
thermal conductivity (Anderson '1998"), e.g., a factor 
of two reduction for FcaSi (Manga ct al. 1996). A low 
value of k is an alternative m echanism to radiogen ic 
heat to avoid core stratification (jStacev fc LopeHl2007f ). 

Initial core temperature: We presume that the cores 
of more massive planets will be initially hotter because 
of the conversion of a greater amount of gravitational 
potential energy into heat. A hypothetical planet with 
a small core that grows concurrently by the descent of 
iron from a cold surface will have a final temperature 
To ~ (3GMp)(10i?pCp) - 22100(Af/Me)2/3 K. These 
core temperatures are unrealistic because much of the 
gravitational potential energy will be dissipated in and 
carried back to the surface by a partially molten silicate 
mantle. Instead, we set the initial temperature at the 
top of the core relative to the melting temperature 
of Mg-perovskite (MgSiOs), effectively the liquidus 
temperature of the mantle, at the CMB. Core tem- 
peratures reaching this would result in the formation 
of a massive magma ocean whose low viscosity would 
quickly cool a hotter core. Alternatively, our initial 
state can be considered the epoch when any such 
magma ocean solidifies. The experimental data and ab 
initio calculations o f the MgSiOa melting temperature 
(jStixrude et al.l [2009^ show that r is well described to 
136 GPa by the relationship r w 5400(P/136 Gpa)"-367. 
This is equ ivalent to a Lindeman law and 7 = 1.44, c.f. 
lAkaogi fc lto (1993), and we simply extrapolate this to 
higher pressures. The mantle is given its steady-state 
temperature as its initial temperature. 



Figure |3] plots the predicted inner core size and the 
surface magnetic field of the Earth at 4560 Myr for 
different values of the initial core temperature (from 
100 K above to 500 K below the MgSiOa mehing 
temperature), critical Rayleigh numbers between 460 
and 3610, and three different values of the thermal 
conductivity (30, 35, and 40 W m-^K "i). We calculate 
a present surface-average terrestrial field strength using 
a dipole moment of 8 x 10^^ A m^. To compare 
with the partially non-dipolar field of the Earth, we 
multiply our model predictions by 80% (Lowes 2007). 
Because 80% is the measured dipole fraction at the 
surface, and the fraction at the core will be smaller, 
our predictions should over-estimate the strength of 
the Earth's current field. At fc > 45 W m~^ K~^ 
the core becomes partly stratified, inc onsistent with 
observations of magnetic field variability (Gubbins 2007t 
[Stanley fc Mohammadi 2008 ). At lower values of fc, the 
present core size and magnetic field strength are can be 
reproduced by combinations of a range of critical CMB 
Rayleigh numbers and initial core temperatures. We 
adopt fc = 35 m~^ K ~^, Ra^, — 1100, and an initial 
core temperature of 275 K below the MgSiOa melting 
temperature at CMB pressure. We use these parameter 
values for all subsequent calculations. 

Figure |4] shows the calculated 4.56 Gyr evolution 
of a 1 M0 planet. The mantle potential temperature 
decreases by 200 K, consistent with paleothermometry 
based on the chemistry of Archea n and Proterozoic 
non-arc basalts (|Herzberg et al.|[2010[) . The temperature 
contrast across the CMB decreases until the onset of 
inner core formation, then rises slightly, as the core 
is heated by latent heat from the inner core, to a 
present value of aroun d 1600 K, only sl ightly higher 
than current estimates (iTateno et al.ll2009| ). Calculated 
mantle heat fiow, radiogenic heat flow, and co re heat 
flow are all consisten t with current estimates ([Buffettl 
l2002HLaveraIll2008D . but the Urey number at 4.56 Gyr 
is ^0.8, signifi cantly higher than in some models 
(jKorenagall2008r ) . The predicted evolution of the surface 
field (or the equivalent d ipole moment) i s comparable 
to that of other models (jLabrossd l2007t lAubert et al.l 
120091: iBreuer et ai]|2009D . It also does not conflict with 
paleomagnetic data, represented by a filtered sample o f 
the lAGA paleointensity database (|Biggin et al.l l2009| ) 
as well as three recent measurements in 3.45 Ga rock s 
(jHerrero-Bervera fc Moizsisll2009HTarduno et al.ll2010( ). 



4.2. Inner cores versus iron snow 

Whether a core grows outward from the center of the 
planet, or inwards from the top of the core depends on 
A, the ratio of the pressure-temperature slopes of the 
adiabat and solidus. This value will vary with depth in 
the core, but is prudently evaluated at the inner core 
boundary or, when no core is present, at the center of 
the planet. For an adiabatic temperature profile. 



A = 



2(37-l)iri(l + ln(p,/po)) 
j[dP/d{p/p,)] 



(46) 



The zero-pressure compressibility Ki appears as a linear 
factor in the EOS, therefore A ultimately depends only 
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on the Griineisen parameter 7, the pressure-derivative 
terms in the EOS, and the pressure at the inner core 
boundary. Because gravity is comparatively low in the 
core, the pressure is relatively insensitive to the size of 
the core, and the total mass of the planet is the most 
important parameter. A will decrease with mass due 
to the increasing incompressibility of Fe (including the 
TFD effect) at high pressures. A < 1 for planets larger 
than 2.5-3 M® (Figure [5]). 

Figure [5] shows the logarithmic ratio of the tempera- 
ture to the melting point versus depth in the cores of 
1-10 M.Q planets. An unimportant constant has been ig- 
nored in each case. The minima in these curves is where 
iron will first solidify. Above 2 M^ the minimum shifts 
from the center (corresponding to inner core growth) to a 
position near or at the CMB. In these cases, iron "snow" 
will freeze out at or near the top of the core. Except 
under an exceptionally narrow range of conditions when 
the temperature coincides with the solidus throughout 
much of the core (Figure |6]) , iron crystals or slabs will 
first form at or near the top of the core and sink, and the 
resulting release of light elements and transport of latent 
heat will stabilize th e liquid core and shut down convec- 
tion (|Williamsll2009D . This effect has been discussed for 
several small bodies in the Solar System (iHauck et al 



20081) 



20i06t iStewart et al|[2007t iChen et ai.|[2008HBland eTa! 



Because there is no net release of latent heat, the 
core will continue to cool and iron will crystallize over 
an expanding depth range (Figure |5]). Double diffu- 
sive convection can drive a "fingering" instability in a 
stratified fluid where two gradients, e.g. a thermal gra- 
dient and compositional gradient, oppose each other, 
but are associated with substantially different diffusiv - 
ities (jHuppert fc Turned llQSTI: iHuppert fc Sparkslll984[ ). 
Double diffusive convection has b een included in a re - 
cent model of Mercury's dynamo (jManglik et al.|[2010[ ). 
In our case, the release and absorption of latent heat 
by crystallizing, sinking, and melting grains of iron will 
create a slightly subadiabatic temperature gradient, and 
buoyancy forces due to temperature will not play a role. 
Double diffusive convection will only occur if (a) the 
diffusion of the light elements is faster than the iron 
crystals, and (b) the crystals are well-coupled to the 
fluid, i.e. there is a high density of nucleation sites and 
crystals are small. Even if convection ensues, it will 
do so in discrete layers each of whi ch has a Rayleigh 
numb er of order the critical value (jHuppert fc Sparksl 
I1984D . Boundary layer theory for finite-amplitude con- 
vection ([Turcotte fc Schubertil200"^ shows that the mag- 
netic Reynolds number of each layer will be of order 
(K/A)Ra^/^, or ^ 10"'^, well below that required to sus- 
tain a dynamo. If the crystals are not coupled to the 
fluid, no convective instability will occur, and instead 
crystals will steadily settle through the fluid. It is not 
clear that a dynamo occurs under these conditions, ei- 
ther. 

4.3. Variation in the dynamo history of Earths and 
super-Earths 

For a given total mass, the mass of a planet's core 
may differ as a result of variation in the ratio of Fe to Si 
in the primordial disk (Neves ct al.,,2009.) . the efficiency 



of core formation (jElkins-Tanton fc Seagerl 2008 ) or im - 
pact stripping of the mantle ( Benz et al.l 119881 I2007D . 
The ability of impacts to remove substantial amounts of 
the mantle from Earth-mass and larger planets might be 
limi ted, and thu s massive Mercurys may not be plausi- 
ble (jMarcus et al. 20 10). Figure [^ shows the predicted 
dynamo history of Earth-mass planets with core sizes 
between 0.5 and 1.5 times that of the Earth. The aver- 
age magnetic field is higher on planets with larger cores, 
mostly due to the smaller r~^ attenuation of the dipole 
field outside the core. 

The initial abundance of heat-producing radioactive 
nuclides (^^K, 232Th, ^ssy, and ^s^U) may also vary be- 
tween planets and between systems as a result of chem- 
ical and thermal processing of solids in planet-fo rming 
disks, and galactic chemical evolution ([Kite et al.| [2009'). 
K is a volatile element and the abundance of ^°K will 
vary with the location of a planet's accretion zone, with 
those planets (e.g. Mars) forming further out having a 
higher abundance (ILodders &: Feglevlll997D . The Galac- 
tic abundances of these isotopes relative to major planet- 
forming elements such as Si are predicted to have de- 
creased substantially with time, and planets of a given 
age may accrete with substantially different inventories, 
but the difference between planets of differe nt ages ob- 
serve d at the present time is much smaller (jKite et al.l 
I2009f ). Figure [71d shows the predicted dynamo evolution 
for planets with 50% and 150% of the Earth's estimated 
initial inventory of radionuclides. The predicted surface 
field is slightly lower with a higher abundance of radionu- 
clides, but is relatively insensitive to this parameter 

Figure [5] plots the predicted dynamo history of 
1-10 M© Earth-like planets with and without plate 
tectonics, and surface temperatures of 288 K (Earth- 
like) and 1500 K (the solidus of basalt). Planets 
much larger than 1.5 M0 with stagnant lids do not 
produce dynamos and even more massive planets have 
lower mantle temperatures that exceed the estimated 
peridotite solidus. Correct modeling of planets with 
internal magma oceans is beyond the scope of this study 
and we terminate those runs. Only planets smaller 
than 2.5 Mq form inner cores (< 2 Mq for those with 
stagnant lids), and the elapsed time before an inner 
core appears, corresponding to a rapid increase in the 
strength of the magnetic field, increases with planet 
mass. The strength of the magnetic field at any given 
age decreases with planet mass. On Earth-mass planets 
with elevated surface temperatures, the strength of the 
magnetic field is higher and the interval before the onset 
of core formation is briefer. 

Planets more massive than 2 M^ (like CoRoT-7b) have 
long-lived dynamos produced by core cooling, but only 
if they have a mobile lid. The duration of the dynamo 
decreases with planet mass and increases with surface 
temperature. The bottom right-hand panel of Figure [8] 
extends these trends out to a mass of 10 M®. 

5. DISCUSSION AND CONCLUSIONS 

5.1. Interpretation of trends and implications 

For an inner core to form, enough sensible heat must be 
removed for the center temperature to reach the solidus. 
We obtain an approximate timescale tc for this by setting 
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Qs — ~Qk and performing the appropriate integration. 
This gives 



of the field at the surface wiU decrease as R^ ^ and thus 



_ AfeA5 f^D^^^ ro(Fe) 
'^ %^kRc V Rl ''rc(pv) 



(47) 



where the sohdii for Fe and perovskite are at central and 
CMB conditions, respectively. For fc = 40 W ni~^ K~^, 
and Re ^ Mp (Figure [2]), the first factor is w 

2Mp Gyr. The second factor is of order unity: Both 
solidii increase with planet mass and so we expect tc to 

increase slightly faster with planet mass than Mp . This 
is born out by our detailed calculations, i.e a 2 M^, planet 
takes more than twice as long to form an inner core than 
an Earth- mass planet (Figure [8]). 

A maximum time ts for the complete solidification of 
the entire core, after the inner core starts growing, is ob- 
tained by dividing the total latent heat by the heat flow. 
(In our cases, the heat flow through the mantle tends 
to be close to or not much larger than the conduction 
through the core itself.) In this case, ts is dependent on 
mass only through the pressure-dependences of certain 
material properties of the core: 



t. 



CpAS 
2kaG' 



(48) 



For standard values, tg ^ 50 Gyr. The iron cores 
of rocky planets with Earth-like composition will never 
completely solidify, and solid inner cores never include 
more than about half of the total core mass in 10 Gyr. 
Thus, planets that do form inner cores (<2 M^) will have 
dynamos that last longer than the main-sequence life of 
the parent star. 

In planets with masses of 2-3 M^, inner cores do 
not grow but instead iron condenses as "snow" at the 
CMB, shutting down convection and the dynamo within 
10 Gyr. In yet more massive planets, the core never 
cools sufficiently for iron to condense. These planets 
have persistent dynamos produced by core cooling. 
However, the temperature contrast between the core 
and the mantle decays, and the heat flow across the 
CMB eventually falls to the core conduction value. At 
that point, (j} = and the dynamo ceases. 

At a given planet mass, the surface field Bp increases 
with increasing core size, and a solid inner core formation 
forms earlier (Figure [7^) ■ The former is a consequence 
of the cubic distance attenuation law and a shallower 
CMB, as well as increased heat flow across a thinner 
mantle. The latter is a result of more rapid cooling, but 
also higher pressure and a higher solidus at the planet's 
center. Dynamo history is much less sensitive to the ini- 
tial abundance of long-lived radionuclides in the planet's 
mantle (Figure [7b) because of the self-reg ulating nature 



of viscosity-dependent mantle convection (jTozei 



ig natu r 
:rl ll972D . 



Bp generally decreases with planet mass. Equations [2] 
and E] show why: The intensity of the field in the core 
Be depends on the available entropy production per 
unit mass (j) and only weakly on core size. All else being 
equal, (p will decrease with mass because the core of a 
larger planet has less surface area per unit mass and a 
thicker mantle through which to lose heat, li cj) ^ R^^, 
this makes Be insensitive to planet mass. The intensity 



will roughly scale as Mj 



^3/4 



Higher surface temperature leads to elevated mantle 
temperatures and lower viscosities. Although the former 
lowers the temperature contrast across the CMB, the 
latter more than compensates because of the exponential 
decrease of viscosity with temperature. The outcome, 
increased magnetic field strength and dynamo lifetime, 
is clearly seen in Figure [H 



5.2. CoRoT-7b 

Our model predicts that CoRo T-7b has a dynamo at 
its estimated age of 1.3-2.2 Gyr (jLeger et al.l 12009:). but 
only if it has a mobile lid. In this case, the dynamo 
is driven entirely by core cooling, and no inner core is 
present. In fact, it is the absence of iron solidification 
(as "snow") that allows a dynamo to persist. Whether 
or not plate tectonic s opera,tes on CoRoT-7b is de- 
bated (iValencia et al.l [20071 : lO'Neih fc Lenardid [20071 : 
iKite et al.ll2009[ ). and the presence of a dynamo may be 
one means of testing this. However, a core dynamo is 
not the only mechanism that can produce a magnetic 
field on such an object (see ^5. 



5.3. Uncertainties and neglected effects 

High-pressure material properties: The properties of 
silicates and Fe at pressures of hundreds or thousands 
of GPa are very uncertain. Although equations of state 
become dominated by electron degeneracy pressure at 
pressures well above 1 TPa, other parameters, namely 
k and a, are very importantjn determining whether 
a dynamo operates. lAittal (|2006l ) has proposed that 
iron has a tricritical point at around 800 GPa at which 
separate solid and liquid phases disappear and the latent 
heat of fusion vanishes. These conditions are reached 
at the center of a 2.5 M^ planet and at the CMB of a 
6 M© planet (Figure [2]). 

Mantle structure: We assume a fully convecting man- 
tle: A layered mantle will sup port a higher temperature 
contrast for the same heat flow (jMcNamara fc van KekenI 
'2000) , and this would tend to suppress dynamo op- 
eration. The lowest 150 km of the Earth's mantle 
(within the D" layer) may b e in the stability field of a 
post-perovskite (ppv) phase (jOganov fc Ono|[2004D and 
ppv may dominate in sufficie ntly cool "super-Earth" 
mantles (jValencia et al.l |2009[ ). If a magma ocean 
existed at the base of the Hadean terrestrial mantle, 
release of latent heat during crystallization would have 
suppressed the coohng of Earth's core (■Labrosse..2007.) . 
Our calculations indicate that the mantles of planets 
with M > 2.5M0 and stagnant lids are partially molten. 
We do not attempt to model the thermal evolution of 
such objects but they are (at least initially) unlikely to 
support a dynamo. 

Other effects of light element in the core: Our model 
does not include melting-point depression in the core by 
light elements, particularly sulfur (S). The effect of these 
is to delay inner core formation relative to the pure Fe 
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case. In that sense our calculations are optimistic in 
predicting dynamos. (This is in contrast to the case 
of small bodies where S can maintain a dynamo by 
delaying the freezing of the entire core). However, we 
choose our initial core temperature based on the present 
size of the Earth's inner core and thus, relative to Earth, 
this calibration may cancel some of the effect of S. We 
include the production of buoyancy flux resulting from 
the exclusion of light elements from a growing solid core. 
However we assume that light elements are uniformly 
mixed into the outer liquid core and do not include 
possible variations in light element concentrations with 
planet mass. The liquid silicate-liquid metal partition 
coeffi cient of sulfur increases with pressure ()Li fc Feil 
I2003D but no experimental data is yet available for 
pressures exceeding those at Earth's CMB. If the trend 
continues to higher mass (and for other elements), the 
power available for dynamos in the cores of super-Earths 
may be smaller than our model predicts. 



Non-radiogenic sources of mantle heat: Dissipation of 
tidal strain could be an additional source of heat in the 
mantles of non-sy nchronously rotati ng planets close to 
their parent stars ([Barnes et al.ll2010l ). The timescale for 
rotational synchronization is a few Myr ([Guillot et al.l 
119961) while that of orbital circularization can be Gyr 
( Rasio et al.lll996l) . but both are strongly dependent on 
semimajor axis. Sufficiently intense heating during this 
period could lead to a long-lived thermal anomaly in a 
mantle that influences the transport of heat much like a 
stagnant lid. We also ignore the p otential complication 



of heating by late giant imp acts ([Roberts et al.l 120091 : 
lArkani-Hamed fc OlsonllMot) . 

Effect of the stellar magnetic field: In the absence of an 
ionosphere, the ambient magnetic field of the host star 
B^, could reach the planet's core: A sufficiently strong 
applied field will inhibit convective motion in transverse 
directions. This will occur when the Stuart number 



St = ^>1. 

pv 

Using the definition of Re^ we re-write this as: 



C7S!i?2 



pX 



> Ren 



(49) 



(50) 



For typical values ct ^ 5 x 10^ S m ^, p ^ W^ kg m, 
A ~ 2 m^ sec~^, L ^ 3 x 10^ m, the condition becomes 

i?* > 60Rem nT. For Re„i ^ 10^ the requirement 
becomes B^ ^ 2 pT. B^ ~ 100 pT near the surface 
of the Sun: Assuming that the stellar magnetic field 
decreases with the square of the distance (because of 
interplanetary plasma), St ^ 1 only within 0.1 AU. 
Thus inhibition of convection could occur in the cores 
of the very closest planets, although the increased 
temperature and conductivity of surface rocks at such 
proximity to the star may block the stellar field (see 



5.4. Other sources of magnetic fields 

In a planet is embedded in an ambient magnetic field 
from its star that varies with orbital period P, eddy cur- 



rents will generate an induced field if a layer is suSiciently 
thick [d) and coi iducting. The co ndition for a substantial 
induced field is ([Stevensonll2003D 



dcr> 



1 



2TTpoRp 



(51) 



For an Earth-size planet where d is in km 
and P in days, the required conductivity a is 
> 2Pd-i(i?p/i?e)-i S m-i. The conductivity of 
anhydrous molten silicates is 0.1-1 S m~^ while that of 
molte n carbonates is three orders of magnitude higher 
(Gaill ard et al.ll2008[) . A planet with a magma ocean 
only a few km deep on a planet with a very short-period 
orbit such as CoRoT-7b (0.85 d) could have an induced 
field. 

Uranus and Neptune have magnetic fields that are 
thought to be generated in a layer of high pressure water 
(ice) (Cavazzoni et al. 1999). The ionic conductivity of 
planetary ices is about 2 x 10^ S m~^ at pressures above 
40 Gpa, condit ions that are reache d atO.7-0.8 radii in 
these planets (IStanlev fc Bloxhaiiil [2006). We added 
liquid water and high-pressure ice (Ice VII) layers to the 
interior model described in H3.1l to calculate the pressure 
at the base of an upper water/ice mantle as a function 
of total planet mass and water mass fraction (Figure [HI . 
A basal pressure of 40 Gpa is reached on planets of a 
few M^ and ice fractions of ~ 10%. Such an envelope 
of water could be easily accomodat ed by the observed 
radiu s (2.7 i?©) of GJ 1214b ([Charbonneau et all 
120091) but not by th e mass and radius of CoRoT-7b 
([Valencia eraIll20To[ ). 



5.5. Conclusions 

Our model of a planet consisting of a silicate mantle 
and pure Fe core reproduces the gross interior structure 
and thermal history of t he Earth. Us i ng the magnetic 
field strength scaling of lAubert et al.l (|2009), the pre- 
dicted evolution of the field strength is also consistent 
with paleomagnetic inferences and the present value. 
When applied to Venus, with a stagnant lid, an elevated 
surface temperature, and enhanced mantle viscosity, 
our model corrctly predicts the absence of a dynamo 
at present. We are thus confident in extrapolating the 
model to rocky planets with approximately Earth-like 
composition but with different mass, surface tempera- 
ture, and/or mode of convection. 

The model predicts that 1-2 Mq planets with plate 
tectonics will sustain a dynamo for up to 7 Gyr by core 
cooling before an inner core forms and the dynamo in- 
tensifies at a time that depends on surface temperature, 
mass, and the relative size of the core. Planets in this 
mass range with stagnant lids will have dynamos only 
once an inner core begins to grow. Planets with higher 
surface temperatures can more easily sustain a dynamo 
because the resulting higher mantle temperatures lead 
to a lower viscosity and hence more rapid mantle 
convection and core cooling. 

Planets larger than about 2 M® do not develop solid 
inner cores because the adiabat is steeper than the 
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solidus. In such planets, iron freezes first at or near 
the core-mantle boundary. The latent heat and any 
light elements released at the top of the core will not be 
able to drive a long-lived dynamo. In planets between 
2-3 M0, this leads to the premature shutdown of any 
dynamo within 10 Gyr. 

In planets with masses > 3 M^, iron never begins 
to freeze in 10 Gyr. If these planets have plate tec- 
tonics, our model predicts that they have relatively 
weak dynamos driven solely by core cooling but which 
can last several Gyr, depending on mass and surface 
temperature. The mantles of planets with stagnant 
lids in this mass range are hot enough to prevent core 
cooling and dynamo operation, but they may also be 
partially molten. This last scenario is not addressed by 
our model. Our model predicts that CoRoT-7b may 
have a dynamo, but only if the planet has a mobile lid, 
and the strength of the field at its surface will be less 
than that at of the present Earth (^ 50/xT) . Advances in 
our knowledge of the properties of planetary materials 
at very high pressure will improve our ability to make 
more quantitative predictions about the intensity and 
longevity of core dynamos in massive rocky cxoplancts. 
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Figure 1. Schematic of our idealized planet. 
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Figure 2. Predicted radius of an Earth-like planet Rp, and its iron core Re, and the central and core-mantle-boundary pressures, Pq and 
Pc, as a function of total mass. 
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Figure 3. Predicted inner core radius and surface field for 4560 Myr-old Earth-mass planets with different values of initial core temperatures 
Tc, critical Rayleigh number Ra* (which parameterizes the uncertain efficiency of heat transfer across the CMB), and thermal conductivity 
of Fe k. Each curve represents a series of model runs with a given Ra* of (top to bottom) 460, 600, 800, 1100, 1560, 2300, and 3610. Each 
point along a curve is represents a run with a different value of Tc- The far right point represents a value 500 K below the lower mantle 
liquidus (5400 K). Each point further to the left represents an increase of Tc in 50 K steps until 100 K above 5400 K. Solid points represent 
runs where the entire core is convecting at 4560 Myr. Open points represent partially stratified cores. All runs assume plate tectonics and 
a surface temperature of 288 K. The present inner core radius and mean surface field of the Earth is plotted as the large circle in each 
pane. The mean was calculated using a dipole moment of 8 X 10^^ A m^. We assume that the magnetic field is a pure dipole, and thus 
our model overpredicts the mean field strength at the surface of the Earth. 
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Figure 4. Predicted evolution of the Earth and its dynamo. Top to bottom: inner core radius relative to outer core; mantle potential 
temperature Q (solid) and temperature contrast across the CMB ATc (dashed); total surface Qm (solid), radiogenic Qr (dashed), mantle 
cooling Qs (dot-dash), and core Qc (dotted) heat flows, and i^Tc-, the power available for the dynamo ( dash-triple-dot); average surface 
magnetic field. In the last pane, 606 measurements from the lAGA p aleointens ity database UBiggin et al.ir2009.) and three recent measure- 
ments in 3.45 Ga rocks l IHe rrcro-Bcrvcra & Mojzsis 2009; Tarduno c t al.ll20"TOl ) are plotted. Only paleointensity measurements older than 
10 Myr, having an error less than 50%, and based on Thellier or Shaw-Thellier techniques with a pRTM check, are shown. 




Figure 5. Predicted ratio A of the temperature-pressure slopes of the adiabat versus the solidus in Earth-like planets with different total 
masses (Earth units). For A > 1, a solid inner core forms. When A < 1, iron condenses at the top of the core. In the latter case the release 
of latent heat stratifies the core and a dynamo is not expected. Solid line: no inner core; dotted line: 50% inner core by mass; dashed line: 
nearly completely solid core. 
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Figure 6. Logarithmic ratio of temperature to solidus with radius in entirely Uquid (Fe) cores in planets of 1 (top) to 10 (bottom) Earth 
masses. Unimportant constants have been ignored. Fe solidification will occur first at the minimum of each curve. Planets with masses 
less than about 2.5 M^ will form solid inner cores. In more massive planets, we predict that iron "snow" will condense near or at the top 
of the core. The concomitant release of latent heat and the absence of buoyancy forces will stratify the core. 
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Figure 7. Predicted history of the average surface magnetic field of Earth-mass planets with plate tectonics and surface temperatures of 
288 K. Dynamos are assumed to produce pure dipole fields, (a) Top: planets with core masses between 0.5 and 1.5x that of Earth, (b) 
Bottom; Planets with an Earth-size core but initial radiogenic element abundances that are 0.5 or 1.5 X the terrestrial case (heavy solid 
line). 
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Figure 8. Predicted history of the average surface magnetic field of planets with masses of 1-fO Mg and the ~4.8 Mgj planet CoRoT- 
7b. Absence of a line indicates that the dynamo is inoperative. Dynamos are assumed to produce pure dipole fields. Heavy solid lines: 
Earth-like planets with plate tectonics and surface temperature of 288 K. Light solid lines: "hot" Earth with plate tectonics and surface 
temperature of 1500 K. Light dashed lines: Venus-like planets with stagnant lid, lOx elevated mantle viscosity, and surface temperature 
of 1500 K. Heavy dashed lines: "balmy" Venus with a surface temperature of 288 K. For the case of CoRoT-7b, the surface temperature 
is set to 1810 K, assuming efficient redistribution of heat (jLeeer et al.ii2009l) . The thick part of the CoRoT-7 curve spans the range of the 
system's estimated age (1.2-2.3 Gyr). Lines terminate when A < 1 and "iron snow" forms at the top of the core, and are absent when 
no dynamo is present or the simulation was discontinued due to the predicted presence of a magma ocean. The behavior of planets more 
massive than 4 M^ is similar to the 4 Mg case. The bottom right-hand panel is the predicted dynamo lifetime as a function of planet 
mass >4 Mg and the two values of the surface temperature. 




Figure 9. Pressure at the base of a water (ice) mantle in solid planets of varying mass and total water mass fraction. Above 40 GPa, 
water can be a sufficiently good ionic (proton) conductor (~ 2 X 10'* S m~^) to produce an internal dynamo analogous to those thought to 
operate in Uranu s and Neptune. The re is also a minimum temperature requirement that is not considered here. The diameter of GJ 1214b, 
~7 Mq, 2.7 -R® l|Charbonneau et al.|[2b09f ). indicates it has a massive volatile envelope that could easily include the required water mantle. 
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Table 1 












EOS Parameters 


Material 


p leg m-3 


Ki [GPa] 


K\ 


Reference 


Ice/water 


1000 


2.2 




4.0 


Sotin et al. (2007) 


HP ice (VII) 


1460 


23.9 




4.2 


Sotin et al. (20071 


Mg-OIivine 


3222 


128 




4.3 


Sotin et al. (2007") 


Perovskite 


4260 


266 




3.9 


Ahrens ('2000) 


Fe(l) 


7019 


154'' 




4.66 


Anderson & Ahrens (1994) 


Fe(.) 


8315 


165 




4.97 


Vocadlo et al. (2003) 


''Adjusted 40% upwards 


from the 


|And. 


arson 


k AhrensI ||1994|) value. 



Table 2 

Thermodynamic Parameters 



Symbol Parameter [units] 



Material 



Value 



Reference 



Cp 


specific heat capacity [J kg -^ K ''-] 


Fe 


850 


Wanf. et al. (2002) 




" " 


MgSiOs-pv 


1250 


Akaogi & Ito (1993) 


k 


thermal conductivity [W m^^K"-"^] 


Fe 


30-100 


Stacev & Loper (2007) 




u 5; 


MgSiOs-pv 


-6 


Goncharov et al. (2009) 


a 


thermal expansivity [K"-"^] 


MgSiOa-pv 


3x 10-5 


Anderson & Duba (1997) 


7 


Griineisen parameter 


Fe 


1.5 


Labrossc (2007) 




u 5; 


MgSiOs-pv 


1.45 


Akaogi & Ito (1993) 


^pIp 


light element density deficit [%] 


Fe 


10 


liAFei (2001) 


AS 


entropy effusion [J kg~^K~^] 


Fe 


118 






" " 


MgSiOs-pv 


130 


Ito et al. (1971) 


A 


magnetic diffusivity [m-^sec"-*^] 


Fe 


2 


Stevenson (2003) 


K 


thermal diffusivity [m-^sec~^] 


MgSiOs-pv 


1 X 10-* 


Turcotte & Schubert (2002) 






Fe 


~ 10-5 


calculated 


(J 


electrical conductivity [S m~^] 
reference viscosity [kg m~-'^sec~-'^] 


Fe 
MgSiOs-pv 


5 X 10^ 

6 X IQis 


Bi et al. (2002) 


V* 


tuning 


TO 


1-bar melting point [K] 


Fe 


1811 


Anderson & Duba (1997) 




liquidus at Earth CMB [K] 


MgSiOs-pv 


5400 


Stixrudc et al. (2009) 



